A LOWER BOUND FOR GARSIA'S ENTROPY 
FOR CERTAIN BERNOULLI CONVOLUTIONS 



KEVIN G. HARE AND NIKITA SIDOROV 

Abstract. Let /3 g (1 , 2) be a Pisot number and let iJ^ denote Garsia's entropy for the Bernoulli 
convolution associated with (3. Garsia, in 1963 showed that H/^ < 1 for any Pisot (3. For the Pisot 
numbers which satisfy a;™ = x™~^ + + ■ ■ ■ + x + 1 (with m > 2) Garsia's entropy has 

been evaluated with high precision by Alexander and Zagier for m — 2 and later by Grabner 
Kirschenhofer and Tichy for m > 3, and it proves to be close to 1 . No other numerical values for 
Hp are known. 

In the present paper we show that Hf3 > 0.81 for all Pisot /3, and improve this lower bound for 
certain ranges of (3. Our method is computational in nature. 



1. Introduction and summary 

Representations of real numbers in non-integer bases were introduced by Renyi fT9\ and first 
studied by Renyi and by Parry [[T6l[T9l . Let /3 be a real number > 1. A (3-expansion of the real 
number x E [0, 1] is an infinite sequence of integers (ai, 02, as, ... ) such that x = J2n>i '^nP~"'- 
The reader is referred to Lothaire, fllSl Chapter 7] for more on these topics. For the purposes of 
this paper, we assume that 1 < P < 2 and E {0, 1}. 

Let pLp denote the Bernoulli convolution parameterized by f3 on I/j := [0, — 1)], i.e., 

{00 
(ai,a2,...)e{0,lf : J]a,/?-^GE 
k=l 

for any Borel set E C where P is the product measure on {0, 1}^ with P(ai = 0) = P(ai = 
1) = 1/2. Since /? < 2, it is obvious that supp (fi/^) = Ip. 

Bernoulli convolutions have been studied for decades (see, e.g., Peres, Schlag and Solomyak 
[fTTll and Solomyak [|22l), but there are still many open problems in this area. The most significant 
property of iip is the fact that it is either absolutely continuous or purely singular (see lessen and 
Wintner fVT\)\ Erdos showed that if /3 is a Pisot number, then it is singular (see [|5|). No other (3 
with this property have been found so far. 

Recall that a number /3 > 1 is called a Pisot number if it is an algebraic integer whose Galois 
conjugates /i 7^ /5 are less than 1 in modulus. Such is the golden ratio r = ^-y^ and, more 
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generally, the multinacci numbers r^, the positive real root satisfying x"^ = + x"^"^ + 

■ ■ ■ + X + 1 with m > 2. The set of Pisot numbers is typically denoted by 5*. It has been proved 
by Salem that S is a closed subset of (1, oo) (see [|20ll ). Moreover, Siegel has proved that the 
smallest Pisot number is the real cubic unit satisfying x^ = x + 1 - see [21]. Amara, [2J, gave a 
complete description of the set of all limit points of the Pisot numbers in (1,2). In particular: 

Theorem 1 (Amara). The limit points of S in (1, 2) are the following: 

(fl = i)l<(f2<^2<^3<X<^3<^^<---<A< V^r+l < ■ ■ ■ < 2 

where 

the minimal polynomial ofifr is ^r{x) = x^"*"^ — 2x'" + x — 1, 
the minimal polynomial ofipr is "^rix) = x^^^ — x*" — ■ ■ ■ — x — 1, 
the minimal polynomial ofx is X{x) = x^ — x'^ — 2x^ + 1. 

A description of the Pisot numbers approaching these limit points was given by Talmoudi 
f23\. Regular Pisot numbers are defined as the Pisot roots of the polynomials in Table ll.ll Pisot 
numbers that are not regular Pisot numbers are called irregular Pisot numbers. For each of these 
limit points (v?r, V'r or x)^ there exists an e, (dependent on the limit point) such that all Pisot 
numbers in an e-neighbourhood of this limit point are these regular Pisot numbers. The Pisot 
root of the defining polynomial approaches the limit point as n tends to infinity. It should be 
noted that these polynomials are not necessarily minimal, and may contain some cyclotomic 
factors. Also, they are only guaranteed to have a Pisot number root for sufficiently large n. 



Limit Points 


Defining polynomials 




<l>^(x)x"±(x'^-x''-^ + l) 

<l>r.(x)x" ± (x"" - X + 1) 

^^(xjx'^i (x'^ + l)(x- 1) 


A 


^/(x)x" ± (x''+^ - 1) 
^^(x)x"±(x''-l)/(x-l) 


X 


X{x)x'^ ± (x"* + X^ — X — 1) 

Afixjx" ± (X^ - x2 + 1) 



Table 1.1: Regular Pisot numbers 



Computationally, Boyd jSl IH has given an algorithm that will find all Pisot numbers in an in- 
terval, where, in the case of limit points, the algorithm can detect the limit points and compensate 
for them. 

Garsia [9J introduced a new notion associated with a Bernoulli convolution. Namely, put 
Dn{P) = |x G //3 : X = YlakP^'' with G {0, 1}| 
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and for x G Dn{(3), 
(1.1) = #<j (ai,...,a„) G {0, 1}" : X = 



A:=l 



Finally, put 



and 



Pn{x) , Pn(a;) 



Hf^ = - y ^log 



2n 

a;eD„(/3) 



i7« = lim 



71 -^OO 



n log /3 

(it was shown in ^ that the limit always exists). The value Hj3 is called Garsia's entropy. 

Obviously, if (3 is transcendental or algebraic but not satisfying an algebraic equation with 
coefficients { — 1, 0, 1}, then all the sums X]fc=i '^kP"'' are distinct, whence p„(x) = 1 for any 
X e Dn{(3),anAHp = log 2/ log /5 > 1. 

However, if P is Pisot, then it was shown in that Hp < 1 - which means in particular that 
(3 does satisfy an equation with coefficients {0, ±1}. Furthermore, Garsia also proved that if 
Hp <l, then /i^ is singular. 

In 1991 Alexander and Zagier in fV\ managed to evaluate Hp for the golden ratio (3 = t with 
an astonishing accuracy. It turned out that Hr is close to 1 - in fact Hr ~ 0.9957. Grabner, 
Kirschenhofer and Tichy [fTO| extended this method to the multinacci numbers; in particular, 
Hr^ 0.9804, Ht-^ 0.9867, etc. They also showed that Hr^ is strictly increasing for m > 3, 
and Ht-^ ^ 1 as m ^ oo exponentially fast. 

The method suggested in [1] has, however, its limitations and apparently cannot be extended 
to non-multinacci Pisot parameters [3. Consequently, no numerical value for Hp is known for 
any non-multinacci Pisot (3 - not even a lower bound. 

The main goal of this paper is to present a universal lower bound for Hp for P a Pisot number 
in (1,2). We prove that Hp > 0.81 for all such /5 (Theorem[9l) and improve this bound for certain 
ranges of P (see discussion in Remark|7]and Proposition [TOl). 

2. The maximal growth exponent 

Denote by £^„(x; P) the set of all 0-1 words of length n which may act as prefixes of /3- 
expansions of x. We first prove a simple characterization of this set: 



Lemma 2. We have 

Sn{x;P) = I (ai, . . . ,a„) G {0, 1}" I < X - f^a^r' < 



k=l 



Proof. Let (oi, . . . , a„) G £n{x; P)\ then the fact that there exists a /9-expansion of x beginning 

Pi 



with this word, implies Y^'^^P ^ < x < a^P ^ + |— j-, the second inequality following 



fromEr+i«;^r'=<g. 
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The converse follows from the fact that if < y < ^zj^ where y = (3'^{x — Y^^=i <^kl3~^), 
then y has a /3-expansion (a„+i, an+2, • • • )• ^ 

The following lemma will play a central role in this paper. 

Lemma 3. Suppose there exists A G (1, 2) such that #£^„(a;; [3) = 0{X^) for all x G //j. Then 

(2.1) H^>\ogf,j. 

Proof. Let (ai, 02, . . . ) be a /3-expansion of x. Denote by p„(ai, . . . , a„) the number of 0-1 
words {a[, . . . , a'^) such that X]fe=i '^kP~'' = Z]fe=i (^'kl^~^- Then, as was shown by Lalley [fT4l 
Theorems 1,2], 



(2.2) ^p,(ai, . . . , a„) ^ 2/3-^^ P-a.e. (ai, as, . . . ) G {0, if . 



Since p„(ai, . . . , a„) < /5) for x = YTk=i <^kl3 ^, we have y/pn{ai, . . . ,a„) < e^A with 

^ 1, which, together with (|2.2I) . implies (12.11) . □ 

Define the maximal growth exponent as follows: 

971/3 '■= sup limsup \/ jj^£n{x; (3). 

It follows from Lemma [3] that 



(2.3) "^^'°i^W, 

Computing 971/3 explicitly for a given Pisot P looks like a difficult problem (unless P is multi- 
nacci - see Section |6l), so our goal is to obtain good upper bounds for 971/3 for various ranges of 
p. To do that, we will need the following simple, but useful, claim. 

Proposition 4. If^Sn+rix; P) < R - i^Sn{x; P) for all n > n^for some > 1 and some r > 2, 
then Tlfs < </R. 

Proof. By induction. 

Let n > uq, and choose kn such that uq + r{kn — 1) <n < nQ + rk^. Then 

i^£n{x-p) < #^„„+,fc„(x;/3) < 2"«i?'=". 

The result follows from 

lim (2"oi?'='Y/" = lim 2"o/"i?'="/" = R^'"- = </R 

n— >oo 71— >oo 

by noticing that ^ — > and ^ ^ ^ as n ^ 00. □ 

Example 5. For the examples in this paper, we give only 4 digits of precision. In fact much 
higher precision was used in the computations (about 50 digits). Let us consider a toy example 
showing how to apply (|2.3I) to P = P^, ^ 1.6737, the largest root of x^ — 2x^ + x^ — x'^ + x — 1 
(which is a Pisot number). 
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Let US first determine (a^ ; /5* ) , dependent upon x . After that we will determine max {x; (3^). 

For ease of notation, we will denote mn{P) = max P). Hence in this case, we are deter- 

mining m2(/9*). Consider the values of x such that x = ^ + ^ + ■ ■ ■ for initial string (ai, 02). 
We see that 

P - - p p^ p^ p^ P P^ 1-1/p' 

This gives us upper and lower bounds for possible initial strings of (ai, 02). 



(ai,a2) 


Lower Bound 


Upper Bound 


(0,0) 


0. 


0.5300 


(0,1) 


0.3570 


0.8870 


(1,0) 


0.5975 


1.1275 


(1,1) 


0.9545 


1.4845 



Table 2.1: Upper and lower bounds for x for initial strings of length 2 of its /9-expansion 

We next partition possible values of x in = [0, 1.4845] based on these upper and lower 
bound. 



Range (approx) 


Possible initial string of expansion 


X e (0., 0.3570) 


(0,0) 


X e (0.3570,0.5300) 


(0,0), (0,1) 


X e (0.5300,0.5975) 


(0,1) 


X e (0.5975,0.8870) 


(0,1), (1,0) 


X G (0.8870,0.9545) 


(1,0) 


X G (0.9545,1.1275) 


(1,0), (1,1) 


X G (1.1275,1.4845) 


(1,1) 



Table 2.2: Initial strings (ai, 02), depending on a; G (0, 1.4875). 



This immediately shows that m2(/3*) = 2. Hence, by induction, #£^„+2(x; /3*) < 2^Sn{x; /?*), 
whence by Proposition!! 3)1^. < v^- By dU), > \ log^^ 2 ^ 0.6729. 

Obviously, this bound is rather crude, and in the rest of the paper we will refine this method 
to obtain better bounds. One thing we need to do is show how one would use this for an entire 
range of P values, instead of just for a specific value. For instance, in the example above, we 
could show that m2{P) = 2 for all /3 > r = ^^/^ ■ In addition, we will want to show how one 
would do this calculation for algebraic P, where we can take advantage of the algebraic nature 
of p. 



6 



KEVIN G. HARE AND NKITA SIDOROV 



3. The algorithm 

Let us consider our toy example of /5 = /3* again. We see that for each initial string (ai, 02), 
we got a lower and upper bound for possible x = aiP~^ + 02/^"^ + ■ ■ ■ . For example, for 
(oi, 02) = (1, 0) these were approximately 0.5975 and 1.1275 respectively. We then used these 
lower and upper bounds to partition into ranges. We next show that if the relative order of 
these lower and upper bounds is not changed, then the partitioning of J/? into ranges can be done 
in exactly the same way. 

Put (ai, . . . , ak)L = El and (ai, . . . , ak)u = + i.e., [(ai, . . . , 0^)^, (ai, . . . 

is the interval of all possible values of x whose /3-expansion starts with (oi, . . . , a^). For exam- 
ple, (1, 0)l = 0.5975 ... and (1, 0)u = 1.1275 . . . This says that if 

(0, 0)l < (0, 1)l < (0, 0)u < (1, 0)l < (0, l)u < (1, 1)l < (1, 0)u < (1, l)c/, 
then we have 



Range 


Possible initial string of /3-expansion of x 


xG((0,0)i,(0,l)z.) 


(0,0) 


((o,i)i,(o,o)f/) 


(0,0), (0,1) 


XG((0,0)t;,(l,0)L) 


(0,1) 







Table 3.1: Upper and lower bounds for x for initial strings of length 2 of its /3-expansion 

as the equivalent table to Table 12.11 For fixed f3, these (ai, 02, . . . , afc)^ and (ai, 02, ... , ak)u are 
called critical points for j3 or simply critical points. 

For each inequality, there are precise values of (3 for where the inequality will hold. For 
example, knowing that (3 > 1, we get that 

(o.ov<(i.o).^^<i^i±^<^ 

So if /5 > r = 1.618 . . . , then (0, 0)u < (1, 0)^. 

This observation means that we need to determine for which values of P we have (oi , a2)L/i/ = 
{a[, a'2)L/u- We will call these values of 13 the transitions points which will affect mn{(3). 

There are some immediate observations we can make that reduces the number of equations to 
be checked. 

• (cti, 02) L = (a'l, 0.2) L and (ai, a2)u = (a'l, ct-'2)u have the same set of solutions. 

• (ai, 02) L = {o-ii 0-2)1/ has no solutions. 

• If ai < a'l and 02 < a'2 then none of 

(ai,a2)L = {a\,a'2)L 

{ai,a2)L = (oi,02)f/ 

(ai,a2)[/ = {a\,a'2)u 

have solutions in Ijj. 
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The first two observations were used when finding all transition points. The last observation was 
made by one of the referees after all of the computations were completed, and hence was not 
used as a means of eliminating equations to check. 

In our length 2 example again, we need to check (after elimination by the three observations 
above), 

{0,0)u = {OA)l, {0,0)u = {1,0)l, (0,0)f/ = (l,l)L (0,l)i = (l,0)i, 
(0,l)i = (l,0)^ (1,0)c/ = (1,1)l (0,l)c/ = (l,0)i, (0,1)^= (l,l)i. 

Solving all of these equations, we see that the only transition points in (1, 2) for length 2 are 
V2^ 1.4142 and r 1.6180. 

So, given that we know niniP*) = 2, and that we have a transition point at r = 1.618 . . . , 
we can say for all /? G (r, 2) that m2{P) = 2. Using a similar method, we can show that for 
f3 e {V2, t) that m2{(3) = 3, and that for f3 G (1, V2) that m2(/3) = 4. 

It is worth noting that these results do not say what happens when /? = or /? = r. The 
transition points will need to be checked separately. 

There is one not so obvious, but important observation that should be made at this point. It is 
possible for an inequality to hold for /5, where /3 is in a disjoint union of intervals. 

For example, we have 

(0,1,1,1,1)l< (1,0,0,0,1)^; 

for P G (1, 0") U (r, 2), where — a"^ — 1 = 0, with a ~ 1.4656. This means that it is possible 
for mn{P) to not be an decreasing function with respect to p. For example 7715(1. 81) = 3, 
m5(1.85) = 4 and 7^5(1.88) = 3. This phenomenon appears to become more common for larger 
values of n. 

4. Numerical computations 

In this section we will talk about the specific computations, and how they were done. The 
process started with length n = 2, and then progressively worked on 77 = 3,4,5,... up to 
77 = 14. We used this process to find the global minimum for all P G (1.6, 2) minus a finite 
set of transition points. The code for for finding transition points, numerical lower bounds, and 
symbolic lower bounds can be found on the homepage of the first author UJJ. 

• For each length in order, find all solutions /3 to 

(ao, ai, . . . , an-l)L/U = {%} • • • ; 0''n-l)L/U 

subject to the conditions mentioned in the previous section. 

• For each of these solutions, check to see if the transition point is a Pisot number. If so, 
we will have to check this transition point using the methods of Section [51 

• Use these transition points to partition (1,2) into subintervals, upon which 7n„(/5) is 
constant. 

• For the midpoint of each of these subintervals, compute 7n„(/9), 

To compute 7n„(/5), we first consider all 0-1 sequences wi,W2, ... of length n. For each 
of these sequences, find their upper and lower bounds, say {ai,a2, . . . } = {wil,wiu,W2l, 
W2U, ■ ■ ■}■ Here the ai are reorder such that ai < ttj+i for all i. We then loop through each 
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interval (a^, aj+i) and check how many of the Wi are valid on this interval. We keep track of the 
interval with the maximal set of valid Wi. 

It should be noted that the number of times we needed to run this algorithm was rather big. At 
level 14, we had slightly more than 300, 000 tests where we needed to find the maximal set. 

These calculations were done in Maple on 22 separate 4 CPU, 2.8 GHz machines each with 
8 Gigs of RAM. These calculations were managed using the Nl Grid Engine. This cluster was 
capable of performing 88 simultaneous computations. 

After this, we looked at all of these subintervals between transition points, and calculated the 
lower bounds for H13 at the endpoints, to find a global minimum. This gives rise to the main 
result of the paper: 

Theorem 6. If f3 > 1.6, and j3 is not a transition point for n < 14, then Hp > 0.81. 

Remark 7. This theorem is weaker than necessary for most values of p. For specific ranges of 
values of P, we actually get a number of stronger results. 

• For most f3 G (1.6, 2.0) have H,3 > 0.82, (99.9 %), and a majority (51.4%) have Hp > 
0.87. Here "most" is a bit misleading. Almost every P has Hp = log 2/ log/9. Of those 
that do not, there is no result that shows they should be evenly distributed, (and they most 
likely are not). So by "most" we mean that for some finite collection of intervals, that 
make up 99.9% of (1.6, 2.0) that all P in this finite collection of intervals have Hp > 0.82. 

• The minimum occurs near ^ 1.8392, (See Figure[2l). 

• For /3 G (1.6, 1.7) we have Hp > 0.87 (Figured, and for p near 2.0 we have Hp > 0.9 
(Figure IH). 



0.94 
0.92 




0.8 
0.78 
0.76 



1.6 1.7 1.8 1.9 2 



Figure 1 . Lower bound for Hp, for Pisot f3 G (1.6, 2.0) and Pisot Transition points 
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Figure 3. Lower bound for Hjs, for Pisot /? e (1.6, 1.7) and Pisot Transition points 
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0.94 - 




Figure 4. Lower bound for if^, for Pisot (3 e (1.98, 2.0) and Pisot Transition points 



5. Calculations for symbolic (3 

In the previous section, we showed for all but a finite number of Pisot numbers /3 in (1.6, 2) 
that > 0.81. To extend the result to all such P in (1, 2), there are still some of Pisot numbers 
that will need to be checked individually. 

These include the finite set of Pisot numbers less that 1.6 (of which there are 12), and the finite 
set of Pisot numbers that are also transition points (of which there are 427). In particular, we get: 



Theorem 8. For all Pisot numbers /? < 1.6 and all Pisot transition points (for n < lA), we have 
Hp > 0.81. 



Combined, this theorem and Theorem [6] yield 



Theorem 9. For any Pisot (3 we have Hp > 0.81. 
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Minimal polynomial oi p 


Pi sot number 


Length 


Lower Bound tor Hp 


X — X — 1 


1.3247 


17 


.88219 


X — X — 1 


1.3803 


16 


.87618 


5 4 3 1 2 1 
X — X — X + X — I 


1.4433 


1 r 

15 


o no r ^ 

.89257 


X"^ — — 1 


1.4656 


15 


.88755 


X — X — X + X — 1 


1.5016 


14 


.90307 


5 3 2 1 

tAj tXj tXj tXj _L 


1.5342 


1 r 

15 


o no T r 

.89315 


x'^ — x^ — x^ ^ x^ — \ 


1.5452 


13 


.90132 


x^ — 2x^ + x^ — x''^ + X — 1 


1.5618 


15 


.90719 


x^ — X** — x^ — 1 


1.5701 


15 


.88883 


X** — x^ — x'' + x^ — 1 


1.5737 


14 


.90326 


nn~^ fY»5 <-v»4 /-v»<^ /-v»^ ry 1 

%Xj tXj lAj lAj lAj lA-/ 


1.5900 


15 


.89908 


ft V -1 

X — X — X + X — 1 


1.5912 


14 


.90023 



Table 5.1: Lower bounds for Garsia's entropy for all Pisot numbers < 1.6 



As a corollary, we obtain a result on small Pisot numbers: 
Proposition 10. All Pisot 13 < 1.7 have Garsia entropy Hj3 > 0.87. 

There are actually a lot of advantages to doing a symbolic check as compared to the numerical 
techniques of the previous section. Some of these include not requiring high precision arithmetic 
and the combining of equivalent strings, both of which has speed and memory advantages. These 
are described in the example below. 

To illustrate the (computer-assisted) proof of Theorem [8l consider as an example (3 = t the 
golden ratio. As before, we wish to find the 

ai a2 ^ ai a2 l/r^ 
h^<x< V — 



r — 1 . Secondly, as ^ = r — 1 and ^ = 2 — r we get 



r r 1 — 1/r 



But now we can find exact symbolic values for these ranges. In particular, we notice that 



(01,02) 


Lower Bound 


Upper Bound 


(0,0) 





r - 1 ^ 0.618 


(0,1) 


2 - r ^ 0.382 


1 


(1,0) 


r - 1 ^ 0.618 


2r-2 ^ 1.236 


(1,1) 


1 


r ^ 1.618 



Table 5.2: Upper and lower bounds for initial strings of length 2 for x = oir ^ + ^ + . . . 



So in particular, it is possible for x to start with both (0, 0) and (1, 0). But if this is the case 
then X = (0, 0, 1, 1, 1, ... ) = (1, 0, 0, 0, ... ) = r — 1. So it is not possible for x to have an 
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infinite number of expansions starting with (0, 0) and an infinite number of expansions starting 
with (1,0). Similar arguments can be used for the other critical point, x = 1. 

So we can discard the critical points and subdivide the possible values of x into the following 
ranges: 



Range 


Possible initial string of the r-expansion 


X G (0,2-r) 


(0,0) 


X e (2-r,r- 1) 


(0,0), (0,1) 


X e (r- 1,1) 


(0,1), (1,0) 


X e (l,2r-2) 


(1,0), (1,1) 


X G (2r-2,r) 


(1,1) 



Table 5.3: Initial string of r-expansion of x, depending on x. 



This immediately shows that m2(r) = 2. Hence, by induction, #£^„+2(x; r) < 2#£^„(x; r), 
whence < V^. By Q, > ^ log^ 2 = 0.7202100. 

The main advantage of this method comes when we have longer strings. In particular, it is 
easy to see that (1, 0, 0) = (0, 1, 1). This allows us to compress information. 





Lower Bound 


Upper Bound 


(0,0,0) 





5 -3r?5i 0.1459 


(0,0,1) 


2r-3 ^ 0.2361 


2 - r ^ 0.3820 


(0,1,0) 


2 - r ^ 0.3820 


4 - 2r ^ 0.7639 


(0,1,1) = (1,0,0) 


r - 1 ^ 0.6180 


1 


(1,0,1) 


3r -4 ^ 0.8541 


2r - 2 ^ 1.2361 


(1,1,0) 


1 


3-r ^ 1.3820 


(1,1,1) 


2r-2 « 1.2361 


r ^ 1.6180 



Table 5.4: Upper and lower bounds for initial string of length 3 for x = air ^ + ^ + a^r ^ + 



This gives that for x G (r - 1, 4 - 2r) we have the initial string of (0, 1, 0), (0, 1, 1), (1, 0, 0), 
and if X G (3r - 4, 1) we have the initial string of (1, 0, 1), (0, 1, 1), (1, 0, 0). 

Our implementation does not maintain a separate entry for (0, 1, 1) and (1, 0, 0), as they are 
equivalent. Instead, the algorithm stores only one of these two strings, and indicates that this has 
weight 2. For the general Pisot /3, this is checked by noticing that (ai, 02, . . . , a„) is equivalent to 

the same word as (&i,62, •••,&«,) if and only if a„x""^H hai = 6„x"~^+6„_ix""^ + - ■ -+61 = 

Cd-ix^~^ + ■ — h Q (mod p(x)) for some q, with p(x) the minimal polynomial for /5, of degree 
d. Given the large amount of overlapping that we see for large lengths, this will have major cost 
savings, both in memory and time. 
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6. The maximal growth exponent for the multinacci family and discussion 

In this section we will compute the maximal growth exponent for the multinacci family and 
compare our lower bound (12.31 ) with the actual values. 

Let, as above, denote the largest root of x™ — x™^^ — ■ ■ ■ — x — 1 (hence r = T2). Define 
the local dimension of the Bernoulli convolution /i^ as follows: 

r / A 1- ^ogfipix -h,x + h) 

dp{x) = hm 

h-^o logh 

(if the limit exists). As was shown in Lalley fT4l, d^(x) = Hp for /i/3-a.e. x G //? for any Pisot 

(3. 

Notice that it is well known that the limit in question exists if it does so along the subsequence 
h = c(3~^ for any fixed c > (see, e.g., Feng [6]). We choose c= {(3 — 1)"^, so 

(6.1) d,{x) = - lim - log, l^x - j-^, ^ + J^)- 

Let P = Tm for some m > 2. 

Lemma 11. Suppose (3 is multinacci, and put 



£,(x) = lim ^#^„(x;/3). 

This limit exists if and only ifdfs^x) exists, and, in this case, 

2 



Proof. Let x = Xlfcli <^kl3~^ and consider (ai, . . . , an), the first n terms of this sequence. We 
see that 



0'n)L — ^^Ctfc/^ ^ 
k=l 

n 00 



k=l k=n+l 
00 00 

k=l k=n+l 
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and 



oo 

-k 



(ai, . . . ,a„)[/ = ^akP ^ + ^ P 



k=\ k=n+l 

oo oo 

k=l k=n+l 

/3" 



x + 



(3-1 



Further, this true, regardless of which representation (oi, 02, . . . ) of x that we take. Hence, if 
(oi, . . . , a„) e £n(x, (3), then for all a^+i, • • • € {0, 1} we have 

n 00 

, 02, . . . , a„)L, (ai, a2, . . . , fln);/) 



fe=l fc=n+l 



— I 3; — -5 — -, a; + 



This in turn implies that 



(6.2) > 2-"#£„(a;; 



Now put 



^n(x;/?) = <j (ai, . . . ,a„) e {0, 1}" | --^ < a; " E«^/^"' ^ |^ [ ■ 



Our next goal is to prove the inequality 



(6.3) \ <2-^4Sn{x-p). 

Let y e {x — 1^^, X + 1^^) have an expansion y — Yl'k=i ^k(^~^- It suffices to show that 
(ai,...,a„) e Sn{x;P). 
By noticing that — <x — y< we get first that 

^ ifc=l fc=l 
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and further than 



X — 


y 


< 


13-1 




-k 


< 


f3-l 




-k 


< 


13-1 


^ - ELi ^kp' 


-k 


< 


l^k=n+l 


^ - ELi ^kP' 


-k 


< 


\-~\OD 

Z^k=n+1 


X - ELi ^ki3' 


-k 


< 


9/3-" 
^/3-l 



' /3- 
(3-1 



Hence (ai, . . . , a„) G En{x; /3) as required. 
Combining (|6.2I) and (16.31) . we obtain 



whence 



1 ~ If /5" 

log/3 2 log^ #^n(a;; /5) < — log^ f^(3 [ x - - — x + 



< log/3 2 - - log;3 #^n(a;; /?). 



(6.4) 

Notice that (16.41) in fact holds for any p. Now we use the fact that P is multinacci. It follows 
from Feng, [[61 Lemma 2.11] that for a multinacci (3 one has ^/jhJx) ~ A/pn(a;') provided 
|x — a;'| < C/3~" for any fixed C > and any x, x' E Dn{P) which are not endpoints of J/j. 
(Here Pn(x) is given by (11.11) .) 

Observe that 

^Sn{x;(3)= Pn{y), 

y&D„{f3y. 
o<y-x<^ 



y&Dn(f3y. 



In view of the Garsia separation lemma (see jSl Lemma 1.51]), each sum runs along a finite set 
whose cardinality is bounded by some constant (depending on /3) for all n. 



Hence ~ y /?) for all x G (0, and together with yield 

the claim of the lemma. 

□ 

Consequently, for a multinacci P, 

2 

(6.5) inf ci^(x) = log^ — , 
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where = {x G (0, : dp{x) exists}. In [6, Theorem 1.5] Feng showed that 



inf dr^{x) 




This immediately gives us the explicit formulae for the maximal growth exponent for the multi- 
nacci family, namely, 



9Jl 



m = 2 

1 

2m+i m > 3. 



In fact, one can easily obtain the values x at which dJip is attained. More precisely, for (3 = t the 
maximum growth is attained at x with the /3-expansion (1000)°°, i.e., at x = (5 + -\/5)/lol] 

For m > 3 the maximal growth point is x with the /5-expansion (10"^)°°. These claims can be 
easily verified via the matrix representation for given in |l6l, and we leave it as an exercise 
for the interested reader. (Recall that the growth exponent for Pn{x) is the same as for ^8n{x] (3) 
for the multinacci case.) 

Finally, since we know the exact values of the maximal growth exponent for this family, we 
can assess how far our estimate (that is, the smallest value of the local dimension) is from the 
actual value of Hp (which is the average value of (i^(x) for jjp-a.e. x). Here is the comparison 
table: 



m 






2 


0.9404 


0.9957 


3 


0.8531 


0.9804 


4 


0.8450 


0.9869 


5 


0.8545 


0.9926 



Table 6.1: Lower bounds and the actual values for H^-^ 



We see that for m > 3 our bounds are far below H^; moreover, our method cannot in principle 
produce a uniform lower bound for all P better than 0.845. However, as a first approximation it 
still looks pretty good. 

Remark 12. We believe (|6.5I) holds for all Pisot /3 G (1, 2). If this were the case, then (12.31) would 
effectively yield a lower bound for the infimum of the local dimension of yU/3. This may prove 
useful, as, similarly to the entropy, no lower bound for c?^ is known for the non-multinacci /3. 
Furthermore, if one could compute the exact value of 9Jt/3, this would yield the exact value of 

inf ^.£7* df3{x). 



This was essentially proved by Pushkarev ifTSl. via multizigzag lattices techniques. 
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7. Acknowledgements and additional comments 

The authors would like to thank the two referees for many useful suggestions. In addition, we 
would like to communicate a question asked to us from one of the referees, that the authors feel 
would make an interesting question for possible future research. 

In Section 6 - besides the multinacci, could you say something on P = {a + 
^/ + 4), with an integer a > 2? (Maybe using results from Komatsu [fT3l .) Or, 
more generally, on numbers P that are root of a polynomial X" — a„_iX"~^ — 
■ • • — aiX — ao, where ctn-i > ■ ■ ■ > c^i > 1? 

We would also like to mention the recent paper by Feng and the second author [7 ], in which 
the average growth exponent for /^-expansions is studied for the Pisot parameters /3. 
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